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Abstract. Searching for and charactering the non-Gaussianity (NG) of a given field 
has been a vital task in many fields of science, because we expect the consequences of 
• different physical processes to carry different statistical properties. Here we propose 

a new general method of extracting non-Gaussian features in a given field, and then 
use simulated cosmic microwave background (CMB) as an example to demonstrate its 
^ , power. In particular, we show its capability of detecting cosmic strings. 

^ I With the cosmological principle as the basic premise, two currently competing 

Q ■ theories for the origin of structure in the universe are inflation [1] and topologi- 

cal defects [2,3]. Although the recent CMB observations seem to have favored the 
former [4], the latter can still coexist with it. In particular, the observational verifi- 
' cation of defects will have certain impact to the grand unified theory, since they are 

an inevitable consequence of the spontaneous symmetry-breaking phase transition 
in the early universe. In addition to the conventional study of the power spectra 
^ \ of cosmological perturbations, another way to distinguish these models is via the 

' search for intrinsic NG — while the standard inflationary models predict Gaussian- 

ity, theories like isocurvature inflation [5] and topological defects [6] generate NG. 
^ , Here we shall propose a new method of extracting the NG from a given field [7], 

^ \ and then apply it to the CMB, which is arguably the cleanest cosmic signals [8]. 

The new method aims to nothing but removing the 'Gaussian' components: the 
mean and the power. Using an n-dimensional field A(x) as an example, the method 
first Fourier transforms A(x) to yield A(k). Then the power spectrum can be 
estimated as Ck = (|A(k)p)fc/l^", where V"^ is the n-dimensional volume of the 
field and k = |k|. Next we define and calculate (Vk with Ck 7^ 0) 



Apfk) 



A(k)-A(o)5(k)lc,-^/x/^ (i: 



where 5(k) is a Dirac Delta and Pk is a given function of k. Finally, Ap(k) is 
transformed back to the real space Ap(x). Now the field Ap has a mean A equal 
to zero and a power spectrum renormalized to P^- For the simplest case = 1, the 
field A is 'whitened' in the Fourier space, and we shall use the superscript 'W to 




denote such whitened fields. In the real space, this means A = ® D + A, where 
(S> denotes a convolution, and D = J dk'^Cl^^e^^'^. Thus A is now decomposed into 
two parts: the 'Gaussian components', A and D, which carry the information in the 
mean and the power spectrum, and the 'NG component', A^, which possesses all 
the remaining information. Therefore, if A is a Gaussian field, then all samples in 
A^ should appear uncorrelated as pure white noise. Otherwise A^ would contain 
all the non-Gaussian features [7] . We note that the above treatment can be easily 
converted to the conventional multipole transform for the GMB, although we shall 
continue to use the Fourier convention, which is appropriate for small CMB fields. 
In this case, we have £ = k and Q = C^. We also notice that the above new 
method is equivalent to the matrix manipulation dp = P^/^C~^/^d, where d = A, 
and P and C are the two-point correlation matrices specified by Pk (with Pq — 0) 
and Cfc respectively. This is similar but different from the Wiener filtering. 

We now test this formalism using simulations. Figure 1 shows six simulated CMB 
components: (a)-(d) (where (d) contains 5 diffuse points) are non-Gaussian, while 
(e) and (f) are Gaussian. They are then hnearly summed to yield A = J2i^{i), 
with RMS ratios ((a)-(f)) 1 : 1 : 10 : 500 : 1000 : 0.2 (Figure 2 left). We then apply 
our method to obtain the extracted NG signal A^ (Figure 2 right). 

In a second test, we simulate a CMB field of (2°)^ (Figure 3(c)): A^ = [AWp] (g) 
Wo + Anoi, where A = Abg (Gaussian background) +A51SW (string- induced CMB 
[9]; Figure 3(b)) -|-Apnt (point source; Figure 3(a)) with RMS ratios 5:1:2, Anoi 
is a 5% noise, and Wj with i = 'p' and 'o' denote the primary and observing beams 
respectively. The whitened field is shown in Figure 3 (d). 




FIGURE 3. Simulated CMB, the extracted NG signal (d), and the power spectra (right). 



With even more tests, the main observation remains the same: in a field A = 
A(G) + A(NG)) regardless how stronger the A(g) is, the NG features of A(ng) can 
always show up in the whitened field A^ as long as Cfc(NG) dominates Ck(G) within 
a certain range of k. In fact, this can be analytically proved [7]. In addition, the 
NG features of uncorrelated NG components do not mix up in the extracted NG 
field A^, even if some of them dominate the others in power. 

Finally we notice that according to equation (1), in principle we can design a 
'window function' to keep the power only on scales where the NG components 
of a field dominate. However, in general we do not know what these scales are and 
thus taking Pfc = 1 is optimal. This may even enable us to find the NG signals of 
unknown physical processes. We acknowledge the support from NSF KDI Grant 
(9872979) and NASA LTSA Grant (NAG5-6552). 
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